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Abstract 

We explore a simple theoretical model to describe the properties of Bose con- 
densed para-excitons in C112O. Taking into account the exciton - phonon interaction 
and introducing a coherent phonon part of the moving condensate, we derive the 
dynamic equations for the exciton - phonon condensate. Within the Bose approxi- 
mation for excitons, we discuss the conditions for the moving inhomogeneous con- 
densate to appear in the crystal. We calculate the condensate wave function and 
energy and a collective excitation spectrum in the semiclassical approximation. The 
stability conditions of the moving condensate are analyzed by use of Landau argu- 
ments, and two critical velocities appear in the theory. Finally, we apply our model 
to describe the recently observed interference between two coherent exciton - phonon 
packets in CU2O. 



PACS numbers: 71.35,+z, 71.35.Lk 

1 Introduction 

Excitons in semiconductor crystals [0 and nanostructures [[| are a very interesting and 
challenging object to search for the process of Bose Einstein condensation (BEC). Nowa- 
days there is a lot of experimental evidence that the optically inactive para-excitons 
in CU2O can form a highly correlated state, or the excitonic Bose Einstein condensate 
Q,@,@]. A moving condensate of para-excitons in a 3D CU2O crystal turns out to be 
spatially inhomogeneous in the direction of motion, and the registered velocities of coher- 
ent exciton packets turn out to be always less, but approximately equal to the longitudinal 
sound speed of the crystal p. 



Analyzing recent experimental ||,||, M and theoretical @-[10] studies of BEC of 



excitons in CU2O, we can conclude that there are essentially two different stages of this 
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process. The first stage is the kinetic one, with the characteristic time scale of 10 ~ 20 ns. 
At this stage, a condensate of long-living para-excitons begins to be formed from a quasi- 
equilibrium degenerate state of excitons (// 7^ 0, T e fj > Ti att ) when the concentration 
and the effective temperature of excitons in a cloud meet the conditions of Bose-Einstein 
Condensation ]lj. Note that we do not discuss here the behavior of ortho-excitons (with 
the lifetime T ort h — 30 ns) and their influence on the para-exciton condensation process. 
For more details about the ortho-excitons in CU2O, ortho-para-exciton conversion, etc. 
see §,1,0. 



The most intriguing feature of the kinetic stage is that formation of the para-exciton 
condensate and the process of momentum transfer to the para-exciton cloud are happening 
simultaneously. It seems that nonequilibrium acoustic phonons (appearing at the final 
stage of exciton cloud cooling) play the key role in the process of momentum transfer. 
Indeed, the theoretical results obtained in the framework of the "phonon wind" model 
fl8fl)||12j| and the experimental observations f3| , [|J , H are the strong arguments in favor 
of this idea. To the authors' knowledge, there are no realistic theoretical models of 
the kinetic stage of para-exciton condensate formation where quantum degeneracy of 
the initial exciton state and possible coherence of nonequilibrium phonons pushing the 
excitons would be taken into account. Indeed, the condensate formation and many other 
processes involving it are essentially nonliner ones. Therefore, the condensate, or, better, 
the macroscopically occupied mode, can be different from n(k = 0) ^> 1, and the language 
of the states in k-space and their occupation numbers n(k) may be not relevant to the 
problem, see [|13] . 



In this study, we will not explore the stage of condensate formation. Instead, we 
investigate the second, quasi-equilibrium stage, in which the condensate has already been 
formed and it moves through a crystal with some constant velocity and characteristic 
shape of the density profile. In theory, the time scale of this "transport" stage, At tr , 
could be determined by the para-exciton lifetime (r para ~ 13 [is [0). In practice, it 
is determined by the characteristic size i of a high-quality single crystal available for 
experiments: 

At tr ~ e/ci ~ 0.5 - 1.5 [MS < r para , 

where q is the longitudinal sound velocity. 

We assume that at the "transport" stage, the temperature of the moving packet (con- 
densed + noncondensed particles) is equal to the lattice temperature, 

T e ff = < T c . 

Then we can consider first the simplest case of T = and disregard the influence of 
all sorts of nonequilibrium phonons (which appear at the stages of exciton formation, 
thermalization || ) on the formed moving condensate. 

Any theory of the exciton BEC in CU2O has to point out some physical mechanism(s) 
by means of which the key experimental facts can be explained. (For example, the con- 
densate moves without friction within a narrow interval of velocities localized near q, and 
the shape of the stable macroscopic wave function of excitons resembles soliton profiles 
0.) Here we explore a simple model of exciton-phonon condensate. In this case, the gen- 
eral structure of the Hamiltonian of the moving exciton packet and the lattice phonons is 
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the following: 

H = H m (ft,j>) - vP ex (#, 4>) + H ph {u, tt) - vP ph (u, tt) + H int (fti>, d,u k ). (1) 

Here tp is the Bose-field operator describing the excitons, u is the field operator of lattice 
displacements, 7r is the momentum density operator canonically conjugate to u, v is 
the exciton packet velocity and, finally, P is the momentum operator. Note that the 
Hamiltonian ([[]) is written in the reference frame moving with the exciton packet, i.e. 
x — ► x — vt and v = const is the packet velocity. 

2 3D Model of Moving Exciton-Phonon Condensate 

To derive the equations of motion of the field operators (and generalize these equations 
to the case of T ^ 0), it is more convenient, however, to start from the Lagrangian. In 
the proposed model, the Lagrangian density has the form 

-E a ft$ - - |(^( x ,t)) 2 (^(x,t)) 2 - |(^(x,t)) 3 ^(x,t)) 3 + 

2 2 

+^(d t u) 2 - ^djUsdjUs -pvd t ud x u + ^-(d x u) 2 - 

-^ t (x,^(x,i)Vu(x,t), (2) 

where m is the exciton "bare"mass (m = m e + ~ 3m e for lsexcitons in CU2O), 
vq is the exciton-exciton interaction constant (u > that corresponds to the repulsive 
interaction between para-excitons in Cu 2 jnj), p is the crystal density, cr is the exciton- 
longitudinal phonon coupling constant, and v = (v, 0, 0). The energy of a free exciton is 
E g +h 2 h 2 /2m. Although the validity of the condition [T5j na B <C 1 (5b is the exciton Bohr 
radius) makes it possible to disregard all the multiple-particle interactions with more than 
two participating particles in H ex , we leave the hard-core repulsion term originated from 
the three-particle interaction in C, i.e. V\ 7^ and < V\ <C uq. For simplicity's sake, we 
take all the interaction terms in C in the local form and disregard the interaction between 
the excitons and transverse phonons. The packet velocity v = v s is one of the parameters 
of the theory, and we will not take into account the excitonic normal component and 
velocity, (T = 0). 

The equations of motion can be easily derived by the standard variational method 
from the following condition: 

SS = 5 J dtdx£(ftfat), ^(x,t), u(x,t)) = 0. 

Indeed, after transforming the Bose-fields tp^ and tp by 

•0(x, t) — > exp(—iE g t/h) exp(zmt>x/^,)-0(x, t), 
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we can write these equations as follows: 



(ihdt + mv 2 /2)4>(x,t) 



h 2 



= ( - —A + i/ W(x, t) + ^i^ t2 ^ 2 (x, t) )</>(x, £) + a Vu(x, t) ^(x, t), (3) 

(^-^A-2^A + ^)u(x,t) =p-V V(^(x,t)). (4) 

We assume that the condensate of excitons exists. This means that the following 
representation of the exciton Bose-field holds: if) = ipo + 5if). Here ipo 7^ is the classical 
part of the field operator if) or, in other words, the condensate wave function, and Sif> is 
the fluctuational part of if>, which describes noncondensed particles. 

One of the important objects in the theory of BEC is the correlation functions of 
Bose-fields. The standard way to calculate them in this model (the excitonic function 
<V ; ( X , t)V^(x', t') >, for example,) can be based on the effective action or the effective 
Hamiltonian approaches ||16|| . Indeed, one can, first, integrate over the phonon variables 
u, get the expression for S e f[(ip, if)^) and, second, use S e fi (or H e s) to derive the equations 
of motion for if> Q , Sif), correlation functions, etc.. 

In this work we do not follow that way; instead, we treat excitons and phonons equally 



||, [17]. This means that the displacement field u can have a nontrivial classical part too, 
i.e. u = u + 5u and u 7^ 0, and the actual moving condensate can be an exciton-phonon 
one, i.e. ^(x, t) ■ u (x, t). Then the equation of motion for the classical parts of the fields 
if) and u can be derived by use of the variational method from £ = C(if),if)*,u), where all 
the fields can be considered as the classical ones. Eventually we have 

(ihd t + mv 2 /2)ip (x,t) = 

h 2 

= ( - — A + i/ |Vb| 2 (x, t) + ^o| 4 (x, t) )Vo(x, t) + a Vu (x, t) ^ (x, t) , (5) 

(d 2 - c 2 A - 2vd t d x + v 2 d 2 x )u (x,t) = p-V V(|^o| 2 (x,t)). (6) 

Notice that deriving these equations we disregarded the interaction between the classical 
(condensate) and the fluctuational (noncondensate) parts of the fields. That is certainly 
a good approximation for T = and T < T c cases [ p~8f| . 

In this article a steady-state of the condensate is the object of the main interest. In the 
co-moving frame of reference, the condensate steady-state is just the stationary solution 
of Eqs. (|H), @ and it can be taken in the form 

V>o(x, t) = exp(-ifd) exp(i<p)0 o (x), w (x, t) = q Q (x), 

where O and q Q are the real number functions, and ip = const is the (macroscopic) 
phase of the condensate wave function. (This phase can be taken zeroth if only a single 
condensate is the subject of interest.) Then, the following equations have to be solved: 

h 2 

A#o(x) =( - —A + z/ o 2 (x) + z^(x) )0 o (x) + a Vq o (x) o (x) , (7) 



4 



( (c 2 - v 2 )d 2 x + c 2 d 2 + c 2 d 2 z )q G (x) = p-VoV^(x). 
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Indeed, the last equation can be solved relative to Vq D . If v < Cj, the corresponding 
solution can be represented as follows: 



Vq (x) = (-y-jj)-^(x) + ^(1-A- 2 ) /^(x-x')-^(xVx', (9) 



where A 2 = cf/(c 2 — v 2 ), and T can be expressed in terms of the Green function of Eq. 
(H). Substituting Vq Q in Eq. (^), we rewrite the latter in the following form: 



where the effective exciton-exciton interaction U e g is induced by the lattice (q Q 7^ 0). It 
can be represented as follows: 



Note that v c g can be positive or negative depending on the value of v. The second term in 
the effective potential (|TT| ) is strongly anisotropic. Moreover, on the cylinder y 2 + z 2 = e 2 , 
its value is negative in the vicinity of x = and positive at the large scales of x (Ox\\ v). 

The possibility of the existence of specific and to some extent unexpected solutions of 
Eq. ( |10|) follows from the fact that the effective two particle interaction between excitons 
can be attractive at small distances between the particles and repulsive at large distances. 
For example, the wave function o (x) may become strongly localized in 3D space because 
of this attraction. 

If v > cj, the formulas for Vq (x) and U e s(yi) are very similar to those obtained in 
the case of v < q. For instance, the solution of Eq. (|8|) can be written as follows: 



(Here A 2 = cf/(v 2 — cf)). Again, the renormalization (|I2|) takes place and a nonisotropic 
part of i7g ff (x) appears in Eq. (|l0|) . However, the effective two-particle interaction between 
excitons remains repulsive, i.e. ^o5(x) + ?7eff( x ) > 0- 



M>(x) =( - — A + z/ ^(x) + J f/ cff (x - x')^(x') dx' + ^(x) )0 o (x), (10) 





Vq (x) = (y " ^ 0o(x) + ^(1 + A 2 ) / T' (x - x'A 2 (x') dx'. (13) 
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3 Effective ID Model for the Condensate Wave Func- 
tion 

Solving Eqs. i n 3D space seems to be a difficult problem (see Eqs. (PH|),(]TT|) ). 

However, these equations can be essentially simplified if we assume that the condensate 
is inhomogeneous along the x-axis only, that is O ( X ) = (fio(x) an d q (x) = (q (x), 0,0). 
Such an effective reduction of dimensionality transforms the difficult integro-differential 
equation (|To|) into a rather simple differential one, and obtained in this way the effective 
f D model for the condensate wave function o °<?o conserves all the important properties of 
the "parent" 3D model. Indeed, iff < q, the following equations stand for the condensate: 

Pxf> {x) = ( -(H 2 /2m)d 2 x + u cS <f)l{x) + U!<f>i{x) )0 o (x), (14) 

d x Qo(x) = const - {a /p(cf - v 2 )) <fil(x), (15) 

where i/ e ff = ^0 ~ &o/ P( c f ~ v 2 )- If ^ > c h Eq. fll~4l) describes the excitonic part of the 
condensate, but with the enhanced effective repulsion u' eS = u + o 2 jp(y 2 — c 2 ). 

The effective two-particle interaction constant v e g is negative if the velocity of the 
condensate lies inside the interval v a < v < q, where 

= \Jcj - (<rg/pi/„) (16) 
can be called the first 'critical' velocity in the model. (Note that outside this interval 

^cff > !.) 

The estimate value of the threshold velocity v Q can be obtained from the following 
formula: 

^ Q \[l - (C c -a) 2 /(87rRy*Ry 7 3 ), 

where C c — C v is the relative volume deformation potential of a semiconductor, (cr — 
C c — C v ), Ry* and Ry are the exciton and atom Rydberg energies, 7 = a B /ai and a\ 
is the lattice constant. (The repulsive exciton-exciton interaction is taken in the form 
Uq ~ 4irh 2 a ex /m). In the case of CU2O oxide, we have v Q ~ (0.5 ~ 0.7)q. 

In this study we will consider the case of v Q < v < q in detail. If the sign of the 
effective two-particle interaction can vary, an extra (repulsive) term should be included 
into the Hamiltonian of a many-particle system to insure the finiteness of a steady-state 
wave function or the absence of wave function collapse in dynamic processes. In our case, 
it is the term z/i(\E'^) 3 ^I' 3 , and V\ > is supposed to be the smallest energy parameter 
in the theory. In the framework of the considered ID model, however, a finite steady- 
state solution of Eqs. (pT4|) , (p~5|) can be obtained without accounting for the "hard core" 
repulsion term V\<$>\. Indeed, we can write out the corresponding solution as follows: 

4> {x) = $ cosh _1 (/3$ a;), d x q (x) = -(cr /p(c 2 - v 2 ))® 2 cosh~ 2 (/3$ x), (17) 

ft = u^ 2 J2 - (a 2 /p(c 2 - v 2 ))$ 2 /2 < 0, (18) 



(3 = (3{v) 



2 cf - V 2 ' 
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The amplitudes of the exciton and phonon parts of the condensate, the characteristic 
width of the condensate, L = (f3(v )<3> ) _1 , and the value of the effective chemical potential 
fx depend on the normalization of the exciton wave function (f) (x). We normalize it in 
3D space assuming that the characteristic width of the packet in the (y, z) plane is finite 
and the cross-section area of the packet can be made equal to the cross-section area S of 
a laser beam. Then we can write this condition as follows: 

|V>o| 2 0, t)dx= S J(p 2 (x) dx = N , (20) 
where N Q is the number of condensed excitons. Immediately, we get the following results: 

*o = L = (^(v))'\ (21) 

,2 _ „,2 . AT . n fc2 



/i = -T^— aUoW) =-^7 L o • (22) 



2 cf - v 2 V25 v '> 2m 

The important "nonlinear" property of the obtained solution (|T7|) is the dependence of 
the amplitudes of exciton and phonon parts of the condensate on the velocity v and the 
number of excitons in the condensate, N Q . For example, 

$1 = & (N , v) ~ N 2 Q (v 2 - v 2 )/(c 2 - v 2 ) 

stands for the exciton amplitude. Notice that the characteristic width of the condensate 
and its velocity are not independent parameters, see (pip . For estimates, the formula for 
L can be rewritten as follows: 

L Z ^ 2a-(na|) 1/2 (^fj'^ 

where n is the average density of excitons in the soliton state. Although fig 1 in this 
formula is multiplied by a small factor, this factor is not small enough to show quantitative 
agreement with the experimentally observed value of L , 2L exp ~ 10 _1 ~ 10 _2 cm. It 
seems to be reasonable that a theory with nonzero temperature (T > will provide a 
more realistic value of the effective size of the packet. 

Returning to the laboratory reference frame, we can write the condensate wave func- 
tion in the form: 

ipo(x,t) ■ u (x,t)5ij = exp ^— i ^E g H — ^ \p,\j t^j exp(imvx)x 

x ^ocosh -1 ^ 1 ^ (Qo - Qotanh^Q 1 ^ - ut))), (23) 

where we count the exciton energy from the bottom of the crystal valence band; 2Q (N , v) 
is the amplitude of the phonon part of condensate and Q Q oc $ D . To calculate the 
energy of the moving condensate within the Lagrangian approach, (see Eq. (0)), we 
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have to integrate the zeroth component of the energy-momentum tensor Tq over the 
spatial coordinates, 

T °(x, t) = E g j>H + W + (W2)(^ W + |(^u) 2 + d jUk d 3 u k + a ^V V u. 

We do not take into account the small correction to this energy due to the quantum 
depletion of the condensate as well as the term (fi/3) 0® in T®. Then the result reads: 

E = Jdr Tq = E cx + E int + E ph = 

( 777/1/ \ / \ ~\~ Tp 1 (J^ 1 

+ (IAI + (-o/3)*=) + 

The value \p\ is a rather small parameter, 

\p\ ~ 6Ry*(na 3 x ) (t; 2 - *; 2 )/(c 2 - V 2 ), 

and the energy of the phonon part of the condensate is estimated as E ph < (5 ~ 6)N \p\. 
However, the back surface of the crystal will experience some pressure when the conden- 
sate reaches this surface and the excitons are destroyed near it. The estimation of the 
maximum value of the pressure in a pulse is as follows: 

p m ~ 10ct $ 2 ~ 1(T 2 ~ 1(T 3 J/cm" 3 

and the main contribution to this value comes from the phonon part. Indeed, one can see 
that the exciton-phonon condensate carries a nonzeroth momentum P ox = P CXiX + P p h,x' 

Pox = J dx{h/2i)(<f)* (x,t) d x (f) (x,t) -d x <f>Z(x,t) <f) {x,t)) - pd t u (x,t) d x u (x,t) = 
= J dxmv $ 2 (x) + pv 

4 Low-Lying Excitations of Exciton-Phonon Con- 
densate 

Although the condensate wave function (j) {x) -q (x) was obtained in the framework of the 
effective ID model, (but normalized in 3D space), we will use this solution as a classical 
part in the 3D field operator decomposition: 

^(x,t) =exp(-ifit)(<f> (x) + 5't(x,t)), (24) 

Wj(x, t) = q (x)5ij + 5uj(x, t), (25) 

where p = p — mv 2 /2. Substituting the field operators of the form fl24|) , (|25|) into the 
Lagrangian density (||), we can write the later in the following form: 

C = £ (e-^V (x), q (x)5 lj ) + £ 2 (^t(x,t),^(x,t), 8u(x,t)) + (26) 



((c 2 -V^ 1{X) ) 
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where C stands for the classical part of £, and C 2 is the bilinear form in the (^-operators. 
As the classical parts of the field operators satisfy the equality 8S \^q,^}q, Uq\ = 0, the 
linear form in the "^-operators" vanishes in fl26|) . 

In the simplest (Bogoliubov) approximation [O,P0|, £ « C Q + C 2 and, hence, the 



bilinear form C 2 defines the equations of motion for the fluctuating parts of the field opera- 
tors. (To derive them we use the variational method: 5S 2 = 5 f£ 2 (5ip^ , ty, 8u) dx.dt = 0.) 
As a result, these equations are linear and can be written as follows: 



n 



2 



iW t <ty(x, t) = (-^A + |/2| + |2z/ - ^ » — I <%(x) + 3M 4 (x) J #(x, t) + 

+ (Ml(x) + 2z/ 1 ^(a;))# t (x, t) + ( x o o (a;)V5u(x, t), (27) 
(<9 t 2 - c 2 A - 2^ + t; 2 (9 2 )5u(x, t) = /9-VoV (0 o (x) W(x, t) + ^(x, t)) ) (28) 



The same approximation can be performed within the Hamiltonian approach. Indeed, 
decomposition of the field operators near their nontrivial classical parts leads to the de- 
composition of the Hamiltonian ([]]) itself, and - as it was done with the Lagrangian - 
only the classical part of H, H Q , and the bilinear form in the fluctuating fields, H 2 , are 
left for examination: 

ff«i?o(V£,Vo, n Q ,u )+H 2 (5ft, 6$, 5tc,5u). (29) 

In this approximation, the Hamiltonian ( p9[) can be diagonalized and rewritten in the 
form: 

H = ^(e-^VoOz), q (x) ) + 8E + J2 a\a s . (30) 

s 

Here, SE is the quantum correction to the energy of the condensate and the index s 
labels the elementary excitations of the system. The operators ar s , a s are the Bose ones, 
and they can be represented by linear combinations of the exciton and displacement field 
operators: 

a s = Jdx ( C/ S (x) <fy(x) + K*(x) <ty f (x) + X sJ (x) 5%(x) + Y aJ (x) ^(x) ) (31) 
a\ = Jdx ( f/;(x) #t (x) + K(x) ^( x) + X * . (x) ^.( x) + y*^) <^(x) ) (32) 



Note that by analogy with the exciton-polariton modes in semiconductors [^T|] the excita- 



tions of the condensate ( ^3|) can be considered as a mixture of the exciton- and phonon- 
type modes, but in this model the phonons come from fluctuations of the uq(x, t)-part of 
the condensate. 

Since the a-operators (see (|30| ) ) evolve in time as simply as 

a s (t) = e~ iu)3t a s a\{t) = e iuJst a\, 

these operators (and the frequencies {oj s }) are the eigenvectors (and, correspondingly, the 



eigenvalues) of the equations of motion (|27j) , (g§) obtained within the Lagrangian method. 
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Then, the time dependent "^-operators" in fl2"T|) , fl2"g|) can be represented by the following 
linear combinations of the a-operators: 



-icjst 



v!(x) ale*"'*. 



6ftx.,i) = ^u s (x) a s e~ 

s 

6ft(x,t) =X)u:(x)dte <w «* + v s (x)a s e- & 

s 

6uj(x, t)=J2 C sd(*) ^~ tulst + CV(x) ate*".*. 



-iuj s t 



(33) 

(34) 
(35) 



Substituting this ansatz (which is a generalization of the u-v Bogoliubov transformation) 
into Eqs. (|27| ) , (|28| ) , we obtain the following coupled eigenvalue equations [[|: 

(L(A) - tiu.) u s (x) + (y Q <fi(x) + 2u 1 <j>t(x)) v s (x) + a o o (a;)VC s (:r) = (36) 

(z/ o 0^(x) + 2u 1 <f>i(x)) u s (x) + (L(A) + fto;,) v s (x) + a o o (x)VC s (x) = (37) 
-p _1 cr o v(0 o (a;) u s (x)j -p~V o v(0 o ( 



x v s x + 



-IUJ, 



vd x 



C,(x) 



0, (38) 



where L(A) = (-^ 2 / 2 ^)A + \p\ + [2u - a%/p{c 2 - v 2 )} (f> 2 Q (x) + 3vi<ft(x). 

To simplify investigation of the characteristic properties of the different possible solu- 
tions of Eqs. (|36"D-(|3"8|), we subdivide the excitations (|33"|)-(|35|) into two major parts, the 
inside-excitations and the outside-ones. The inside-excitations are localized merely inside 
the packet area, i.e. |x| < Lq and (fl{x) ~ const, whereas the ontside-excitations propa- 
gate merely in the outside area, i.e. |x| > (1 ~ 2)L and 4> 2 {x) ~ <& 2 exp(— 2\x\/L ) — > 0. 



4.1 Outside-Excitations 

For the outside collective excitations, the asymptotics of the low-lying energy spectrum 
can be found easily. Indeed, if we assume that 4>„(x) ~ in the outside packet area, the 



equations (|27| ) and fl28|) begin to be uncoupled. Then, Eq. (|27|) describes the excitonic 
branch of the outside-excitations with the following dispertion low in the co-moving frame: 



^ ex (k) w |£| + (Ti 2 /2m)k\ (u k (x)~e Jkx , v k (x) « 0), (39) 

and Eq. (p8|) describes the phonon branch and yields the spectrum uj p h(k) = q|k| in 
the laboratory frame of reference. Then, the exciton field operator, which describes the 
exciton condensate with the one outside excitation, has the following form: 

?/>(x, t) ~ exp(— i(E g + mv 2 /2 — \p,\)t ) exp(imra) o (x — vt) + 

+ exp(—i(E g + mv 2 /2 — |/2|)t ) exp(imt>x) |exp(— + fik 2 /2m + k x v)t ) exp(ikx) u k | 

It is easy to see that such a collective excitation, c<j ex = \jl\ + hk. 2 /2m + k x v, can be 
interpreted as an exciton with the energy E g + h 2 k 2 /2m, where hkj = hkj +mv6ij. Then 
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Figure 1: Moving exciton-phonon condensate, 4> (x — vt) ■ u Q (x — vt)dij, and inside- and 
outside-excitations of the condensate. (Longitudinal exciton-phonon excitations, k 1 1 Ox, 
are schematically depicted.) 

we can compare the condensate energy E (N ) and the energy of the condensate with one 
outside excitation, 

E (N - 1) + E cxc (~k) + E ph (k') « E (N ) - N E o (N o ) + (E g + h 2 ~k 2 /2m) + h Cl \k'\ « 

2 2 2 

« E (N ) + (h 2 ~k 2 /2m - mv 2 /2) + 3(|/2| + (^o/3)^ ) - $g + ftq|A/| 

c 2 11 P\ c i v ) 

> E (N ) in k -> (fc' ^ 0) limit. 

Note that the energy (and the momentum) of the phonon part of the condensate changes 
after exciton emission. We assume that the transformation N —> N — I (or emission of 
an outside exciton) corresponds to the situation when the outside exciton and the outside 
phonon(s) appear together, and the phonon is emitted with the energy compensating the 
changement of — 5E p h in the phonon part of the condensate. 

However, the condensate collective excitations are uncoupled only in the k — > limit, 
i.e. A = 2n/k ^> L . It follows from the structure of Eqs. (|36D-(|38|) that the coupling 
between the outside phonon and exciton branches is originated from the condensate "sur- 
face" area, i.e. from the scale L Q < \x\ < 3L . Indeed, (f) (x) and (f>' (x) ~ ±4> (x)Lq 1 
cannot be put equal zero in this area, and the "particle"- and the "hole" -type components 
of the exciton operators, namely u s ~ e tkx and v* ~ e -ikx^ s ] lou i ( j De both different from 
zero and spatially modulated in the surface area, at least for the excitations with A < 2L . 
We left this question for future investigations. 
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4.2 Inside-Excitations 



To simplify the calculation of inside-excitation spectrum (see Eqs. (^)-(^) ) we will 
use the semiclassical approximation p0 |. In this approximation, the excitations can be 
labeled by the wave vector k in the co-moving frame, and the following representation 
holds: 

u s (x) = u k (x)e^W, v s (x) = v k (x)e^W, C s ,,(x) = C kj (x)e^ x \ (40) 

where the phase ^k(x) ~ <^ +kx, and u k (x), v k (x), and Cj ik (x) are assumed to be smooth 
functions of x in the inside condensate area. Notice that the k- and x-representations 
are mixed here. This means that the operator nature of the fluctuating fields is factually 
dismissed within the semiclassical approximation, However, the orthogonality relations 
between u s and v s i, and, hence, between u k and v k / come from the Bose commutation 
relations between the operators a s and a\, [fL9|l , pPf . For example, Eq. (|33|) is modified as 
follows: 

Ah- 

/7r ^u k (x)e^ k(x) e- ia;k(x) * + v*(x)e- i¥,k(x) e^ k(x)i , (41) 

and the ms«G?e-excitation part of the elementary excitation term in (|30"|) , 
Es • • • « E a , out + Es, surf + E fl , in ■ ■ • , can be written in the form 

V hu s &l& s ~ [ d ^ d *. ^ k (x) n k (x). (42) 

S ,in J ( 27F ) 

Note that the semiclassical energy Tiuj^x) of the inside-excitation mode is supposed to 
be a smooth function of x, (i.e. at least as smooth as <pl(x) ~ const in the "inside" 
approximation) . 

Although the low-lying excitations cannot be properly described within the semiclas- 
sical approximation, we apply it here to calculate the low energy asymptotics of the 
spectrum. In fact, all the important properties of these excitations can be understood 
within this approach. 

Substituting ( fiOD into Eqs. (|36|)-(|38|), we transform these differential equations into 
the algebraic ones ( L(— k 2 ) = L(A — > — k 2 ) ): 



(L(-k 2 ) - hu*) u k (x) + (z/ o 0oW + 2^i0 o (x) j v k (x) + O r o o (x)zkC k (x) = 0, (43) 

(z/ o 0o(x) + 2i/irf(x)) u k (x) + (L(-k 2 ) + hu k ) v k (x) + <r </> o (x)zkC k (x) = (44) 

p" l (TQ(j) {x)ikj u k (x) + p^ao^o^ikj v k (x) + (w k + vk x f - c 2 k 2 C k j(x) = 0. (45) 

After some straightforward algebra, we can write out the equation for the exciton-phonon 
excitation spectrum: 

( (cu k + vk x f - c 2 k 2 ) x 
x \hco k ) 2 - (L(-k) - U^(x) + 2Mi{x))) (H-k) + + 2^(x))' 
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L(-k) - (vo^(x) + 2^(x)) ) 2 ^l(x) (c 2 k 2 ). (46) 

Note that within the semiclassical description of the msitie-excitations, the low energy 
limit means k — > k where ko is the momentum cut-off, 

(h 2 /2m)k 2 ~ \fl\ = (h 2 /2m)L 2 . 

The inequality k > k ~ Lq 1 ensures the function Tia^x) in (fIBl) to be real and positive. 
Indeed, only the excitations with the wave lengths A < (2 ~ 3)L can be considered as 
the inside ones. The presence of the "hard core" terms, const v^^x), leads to a slight 
renormalization of the value of the momentum cut-off. However, this renormalization does 
not factually change the characteristic properties of the possible solutions of Eq. ([461). 
We will mark the "hard core" terms by e + > 0. For example, in the low energy limit, 

k « k + 5k, 5k — > 0, 

the exciton part of the l.h.s. of Eq. fl46|) - i.e. the formula inside the square brackets - 
can be reduced to the form: 

- ~ k ° } + F(x) + e + ) ^ 2(k 2 2 ~ k ° } + F(x) + 2^ 2 {x) + e+ ) , (47) 

where F(x) = [all pic 2 — v 2 ) — u ) ($ 2 — 4> 2 (x) ), and the following two estimates hold: 
F(x) ~ 2\p\{2x/Lq) 2 at x ~ and F{x) ~ 2\p\ at x > ±2L . 

There are two different types of the inside-excitations, the longitudinal excitations 
and the transverse ones. The later have the wave vectors k perpendicular to the x (in- 
direction. Although Eq. ([46"D can be solved exactly for the transverse excitation spectrum 
T2J , taking into account the coupling term changes the values of excitation energies 



slightly, and the excitations can be approximately considered as of the pure exciton or 
phonon types. Then we have the acoustic phonon dispersion low for the phonon branch 
and the following spectrum for the exciton branch: 




(h^ k± y ~ — (k 2 - k 2 ) + e + — (k 2 - k 2 ) + 2vcg{x) + e 



~ 2^ (fc2 ~ fc ° 2) 2V ^ X) + 2u ^ e+ < ( l^i + ( h2 / 2m ) k l ) 2 • (48) 

The smooth function u; fcx (x) > has a gap when k — > k , but unlike the case of the outside 
excitations, the gap value is determined by the "hard core" repulsion term and is much 
less than \fi\. Furthermore, if we let (formally) the x coordinate in ( PD change in the area 
of \x\ > Lq, the dispersion low fwj eKt f. ± {x) reproduces the outside excitation asymptotics, 

+ (h 2 /2m)k 2 _. However, inside the condensate, we obtain a strong deviation of the 
collective excitation spectrum from the simple excitonic one. 

In the case of the longitudinal excitations, the mode interaction is non-negligible in the 
low energy limit. For instance, the "exciton" root of Eq. (EH) with the nonzeroth r.h.s. 
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exists if k x > 2.5 Lq 1 . Therefore there are no distinct exciton and phonon modes, and the 
cases k x > and k x < are different because of different position of the "bare" phonon 
root on the energy axis. Then, on the energy axis %u, the modified phonon spectrum is 
located higher than the phonon frequencies and the modified exciton spectrum is located 
lower then %u^ k . Yet, like the case of transverse excitations, the same inequality and 
the same (formal) asymptotics are valid for the lower branch of the spectrum: 

< u kx < \p,\ + (h 2 /2m)k 2 x . 

The approximate formulas for the longitudinal spectrum are too cumbersome to be pre- 
sented here. However, the phonon component of the excitonic longitudinal excitation, 
Ck,i{x) 7^ 0, can be found approximately from Eq. (|45"D by use of the Bogoliubov form 
for the wave functions u^(x) and v^(x) fl2"0f : 



. , (I \ L[—k 2 ) + huk 2/ ( 1 \ L(—k 2 ) — Tioj k 



where L(—k 2 ) « h 2 (k 2 — k 2 ,)/2m + ucj) 2 (x) + e+. The effective condensate volume V c g ~ 
4SX is used to normalize the u- and v-wave functions of the inside excitations, J dr(\u k \ 2 — 
|vfe| 2 ) = 1. Subsiquently, we get 

r , , , p~ l o Q (p {x)ik x (u k (x) +v k (x)) 

CtMx)ai — ^-(^ +.*,)> — • (50) 



One can use the zeroth approximation, Eq. (f|8j), for u kx in (^9|),(p0|). Then we can 
roughly estimate the maximum value of |Cfc(x)| 2 in the low energy limit (k — > 2.5L ( ^ 1 , 
a; <C c|fea.|): 

|C fc | 2 ~ (\jl\/p(c 2 - V 2 )V eS ) Lq «< Lq- (51) 

To investigate the stability of the moving condensate in relation to the creation of 
inside excitations, we can calculate the energy of the condensate with the one inside 
excitation described by the following set: k, u k , u k and Vfc, and C k . Although such 
an excitation was defined in the co-moving frame, calculations should be done in the 
laboratory frame. Returning to the lab frame, we represent the exciton and phonon field 
functions as follows: 

4> (x — vt, t) — > <p (x — vt, t) + exp(— i(E g + mv 2 /2 — |/2|)t ) exp(imvx) <5\I/(x, t), (52) 

<5*(x, t) = u k (x - vt)e ikx e- i ( Uk+kaV)t + v k (x- vt ) e - l]iX e +i{uJk+k * v)t , 

u Q (x — vt) — > u Q (x — vt) + C k {x — vt) exp(ikx) exp(— i(ou k + k x v)t ) + c.c. . (53) 

As the inside excitation is considered as an fluctuation, the number of particles in the con- 
densate and its energy can be estimated as N Q —J cLx Sip^Sip and E (N )—dNE J dx. Sijj^Sij), 
respectively. 

The zeroth component of the energy-momentum tensor can be represented in the form 
T ° = T O °(0 O , u Q ) + T° m (5&, <ftf , 5u | fa, u Q ), 
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where the first part gives the condensate energy E Q and the second part will give the 

7-0(2) 
-0 ) 



energy of the inside excitation, E- m . After substitution of (|52D,(|5"3"D into E- m = J dxT^ 2 ' 



it can be rewritten as follows: 



2 

x I • 



En = J dx(E g + mv 2 /2 - \fi\)6^6ip + J dxh(uj k + k x v)(\u k \ 2 - |v fe | 2 ) + 2p\C k \ 2 {uj k + vk 

(54) 

The first term appearing in (|54D , 5ip^5ip — *> \u k \ 2 + \v k \ 2 , vanishes in the following formula 
for the total energy: 



E + E in « E {N ) + (2\fi\ + v <S> 2 o ) / dx8ftdi> + J d*h{u k + k x v){\u k \ 2 - 



Vfe 



1 2> 



cf + V 



2 

., , -p-^rQl I dxSftSj + / d*2p\C k \ 2 {u k + vk x f. (55) 
c « + v P\ c i — v ) J J 

Here the last two terms compensate each other approximately, and all the interesting 
effects come form the exciton part of (|55|). Although u k + k x v can be negative, its negative 
value can be compensated by the term 2\p,\ + z/ $q if the velocity of the condensate is close 
to q or the exciton concentration is high enough. Therefore, there is a second critical 
velocity v c in the theory. If the velocity of the condensate is less then the velocity v c and 
v a < v c < Q, the condensate is unstable in relation to creation of the inside excitations, 
i.e. E a + E m (k x ,v) < E D in the lab frame. To estimate the value of v c , we solve the 
following equation: 

(3 ~ 5)hvL^(v) ~ (a 2 /(cf - v 2 )p) <S>l(v), 

which can be reduced to (3 ~ 5)hv ~ (o"o/(cf — v 2 )p) (N Q /2S). For example, for the packets 
with the exciton concentration of n ~ 10 14 ~ 10 15 cm -3 , the estimate is as follows: 

c t -v c a 2 N 0.1 
c\ pcf 2 b nci 

5 Interference between Two Moving Packets 

There are at least two interesting problems that can be examined in the framework of the 
proposed model (Eqs. flUD). The first problem is the investigation of the condensate 
steady-state and its stability, calculation of the low- lying excitation spectrum, etc.. (Some 
part of this program was presented in the previous sections.) The second one is the 
investigation of interference between two moving packets p3| . In this case the problem is 
essentially nonstationary. Indeed, the amplitude and the shape of the resultant moving 
packet are expected to change in time f24|. These effects, however, can be considered 
theoretically by a numerical solving of Eqs. ([5|),(|6]) with the proper initial conditions. 
(We will disregard the influence of noncondensed particles, the condensate depletion and 
nonequilibrium phonons on the dynamic processes being considered.) For example, if two 
"input" packets have the same velocity (vi = v 2 = v) and shape, we can write the initial 
conditions in the form 

V> (x, t = 0) • u (x, t = 0) = o (x) • q G (x) + e^0 o (x + vr) • q Q (x + vr), (57) 
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where 5(p = const is the macroscopical phase difference between the two "input" conden- 
sates, and r = const is the time delay between them. 

In the simplest (quasi- ID) model, the following equations govern the dynamics of the 
two "input" packets: 

h 2 

(ihd t + mv 2 /2)ip (x,t) = (-— ^ + z/ |V ; o| 2 + ^i|V'o| ) ipo(x,t) + a d x u {x,t) ip (x,t) (58) 

2m 

(d 2 - (c 2 - v 2 )d 2 x - 2vd t d x )u (x,t) = p- l a d x \ip \ 2 (x,t). (59) 

Then the initial conditions ( |5"7| ) can be written in the explicit ID form by using the exact 
solution ([H]) of the model (|T^,[I5]). Note that the amplitudes of 4> (x) and du Q (x) are 
defined from the normalization condition and depend on the values of v and N Q , and, 
hence, the amplitudes of the "input" condensates in fl5"7| ) have the same values. 

Some predictions of the form of the expected solution can be made easily. Indeed, 
the shape and other characteristics of the steady-state solution of (0), (0) will depend 
mainly on the value of the parameter x /L where Xq = vt and L is the characteristic 
condensate width. If Xq/L < 1 ~ 2, the nonlinear interaction between the packets 
plays an important role in the process of total wave function formation. (Notice that the 
lower limit of xq, Xq = r*v, is defined by the time scale of formation of a condensate 
wave function.) It is reasonable to assume that in the limit of strong interaction between 
condensates, the N a + N Q excitons can form the single condensate wave function flTTD in 
the steady-state regime. Then this wave function can be written (in the laboratory frame) 
as follows: 

ipo(x, t) ■ u (x, t) ~ exp(—i(p,(2N , v) + mv 2 /2)t ) exp(imvx) x 

x exp(i(p)(j) (x — vt] 2N Q ) ■ q Q (x — vt; 2N Q ). 
As the dynamic equations (|58|),(j59|) conserve the energy, 

E in (N,v; x /L ) = E out (2N,v), 

v cannot be equal to v in theory. Moreover, if the parameter Xq/L is small enough, it 
can be only the approximate equality, v ~ v. 

In the case of xq/Lq 3> 1, one can disregard the influence of the mutual nonlinear 
interaction on the dynamics of the packets. In this approximation, the packet moving in 
the crystal can be modeled by the following formula: 

if) (x,t) ■ u {x,t) ~ exp(-z(/2(iV ) - mv 2 /2)t)(p (x; N Q ) ■ q (x; N a ) + 

+ exp(i5(p) exp(-i(jl(N ) - mv 2 /2)t )(j) (x + x ; N Q ) ■ q Q (x + x ; N Q ). 

An interesting and noninvestigated case in the interference problem is the condensate 
dynamics after posing nonsymmetric initial conditions. In fact, the amplitude and the 
velocity of the "input" packets can be different, for example, N 2 > Ni and t>2 > v±, v 2 ||vi. 
We use here the experimental result |5|] that at v > Vq, the velocity of the condensate 
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depends on the laser power or, equivalently, on the initial concentration of excitons, i.e. 
v = v(N). (Note that in theory the exciton amplitude $ D is the function of N Q and v.) 

If the exciton concentration in the first packet, ni, is close to the Bose condensation 
threshold value and the exciton concentration in the second packet, n 2 , can be made 
^> ni, the velocity difference between condensates can reach (0.2 ~ 0.3)q. Then, in the 
reference frame moving with the first packet, the initial conditions can be taken as the 
following: 

</>o(x, t = 0) • Uo(x, t = 0) = o (x; Nj) ■ q G (x; Ni) + 

+ exp(i5ip) exp(zm5vx) o (x + x ; N 2 ) ■ q Q (x + x ; N 2 ), (60) 

where 5v = v 2 — v 1; x Q = Ujt, and the second packet moves in this frame of reference. In 
the case of such the initial conditions, the regime of strong nonlinear interaction between 
the condensates is unavoidable. Following the logic of the soliton theory, we speculate 
that the steady-state solution may consist of two packets moving with different velocities 
and with different exciton concentrations. Roughly speaking, the "input" packets could 
exchange their places, i.e. the both packets survive after collision, and the first packet 
arrives at the "detecting" boundary of a crystal after the second one: 

^o(x, t) ■ u Q (x, t) ~ exp(-i(/i(iVi) - mvf/2)t ) <p {.x] N x ) ■ q (x; N x ) + 

+ exp(— i(/2(i\T 2 ) —mv\/2 + mSv 2 /2)t ) exp(z5<^) exp(im 5v x) x 

x0 o (x + x - 5vt; N 2 ) ■ q Q (x + x Q - 5vt; N 2 ). 

However, the hypothesis about the solitonic character of packet collisions in 3D needs 
both numerical and experimental evidence. 



6 Conclusion 

In this study, we considered a model within which the inhomogeneous excitonic condensate 
with a nonzero momentum can be investigated. The important physics we include in 
our model is the exciton-phonon interaction and the appearance of a coherent part of 
the crystal displacement field that makes the moving condensate of the exciton-phonon 
one. Then, the condensate wave function and its energy can be calculated exactly in the 
simplest quasi-lD model. We believe that the transport and other unusual properties 
of the coherent para-exciton packets in C112O can be described in the framework of the 
proposed model properly generalized to meet more realistic conditions. 

As the exciton-phonon interaction is very important in any processes involving excitons 



in lattices |26| , we can speculate about the possibility of use of piezoelectrical transducers 
to pump acoustic waves into the system condensate + lattice. Moreover, the transducers 
could be used to register the phonon part of the coherent packet in experiments in which 
the condensate is formed by optically inactive excitons and phonons. 

We showed that there are two critical velocities in the theory, Vq and v c . The first 
one, vq, comes from the renormalization of two particle exciton-exciton interaction due 
to phonons, and the inhomogeneous soliton state can be formed if v > Vq. The second 
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one, v c , comes from use of Landau arguments jllj for investigation of the dynamical sta- 
bility / instability of the moving condensate. Within the semiclassical approximation for 
the condensate excitations, we found the condensate is unstable if v < v c . It is interesting 
to discuss the possibility of observation of such an instability when the condensate can 
be formed in the inhomogeneous state with v ^ 0, but with v < v < v c (n,v). Such a 
coherent packet has to disappear during its move through a single pure crystal used for 
experiments. As the shape of the moving packet depends on time, the form of the regis- 
tered signal may depend on the crystal length changing from the solitonic to the standard 
diffusion density profile. 

We did not concentrate on detailed investigation of excited states of the moving 
exciton-phonon condensate in this study. First, the possibility of their observation is 
an unclear question itself. Second, the stability conditions of the moving condensate - in 
relation to the creation of condensate excitations - can be derived from the low energy 
asymptotics of the excitation spectra at T -C T c . However, the stability problem is not 



without difficulties P2"|,||25||. One can easily imagine the situation when the condensate 
moves in a very high quality crystal, but with some impurity region carefully prepared in 
the middle of the sample. Then the impurities could bound the noncondensed excitons, 
which always accompany the condensate, and could mediate, for instance, the emission 
of the outside excitations. The last process may lead to depletion of the condensate and, 
perhaps, some other observable effects, such as damping, bound exciton PL, etc.. On the 
other hand, the inside excitations could manifest themselves at T 7^ by the effective en- 
largement of the packet length, L — * L eff , T 7^ 0, or by interaction with external acoustic 
waves. 
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